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Radiation  Pattern  Sidelobes  and  Null 
Filling  Produced  by  Aircraft  Vibrations 


1.  INTRODUCTION 

Presently  there  is  considerable  activity  devoted  to  the  development  of  very 
low  sidelobe  antennas  and  antennas  with  movable  nulls.  Such  systems  clearly  can 
work  well  if  ground-based  but  when  mounted  on  aircraft  there  must  be  some  con- 
cern about  the  effect  of  aircraft  vibrations  on  the  radiation  pattern.  In  this  report 
we  have  developed  a general  formalism  to  consider  this  problem  by  considering 
the  radiation  pattern  of  an  arbitrary  array  of  radiating  elements  which  is  subjected 
to  vibrations.  In  Section  3 and  Appendix  A,  we  treat  the  case  of  small  (compared 
with  a wavelength)  amplitude  harmonic  vibrations,  and  in  Section  4 we  study  the 
situation  when  the  vibration  amplitudes  are  comparable  with  or  greater  than  a 
wavelength. 

2.  ANALYTICAL  PRELIMINARIES 

The  scalar  field,  u(r,  t),  produced  at  a point  £ by  a source  distribution  f(s,  t) 
can  be  written*  as 

(Received  for  publication  23  Januarj  1979) 

1.  Panofsky,  W.,  and  Phillips,  M.  (1955)  Classical  Electricity  and  Magnetism. 
Addison-Wesley,  Reading. 
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IUJI1JL  L IUJPM 


«(£,  t)  = ji 


x rffikuM] 

**  JJJ  |r-s| 


d3s  , 


(1) 


where  c is  the  speed  of  light  and  the  integral  is  over  all  space.  We  now  assume 
that  the  source  is  an  arbitrary  collection  of  + Ng  + 1 moving  monochromatic 
point  sources.  We  can  then  write 


Nr 


iu_t 


f(s,  t)  = Re  £ ane  ° 6[s  - £n(t)l  . 


(2) 


n=-N, 


th 


where  R denotes  " real  part  of, " aR  is  the  complex  amplitude  of  the  n point 
source,  s (t)  is  its  position  and  6 (...)  is  the  Dirac  delta  function.  If  we  assume 


that  the  source  motion  is  non-relativistic,  we  may  approximate 

i(t)1  * 

Then,  upon  inserting  (2)  into  (1)  we  get 


{i-fnO  ' ~c“)]  -ala-an( 


N„ 


u(r,  t)  a Re  X an 


n = -N, 


exp  jiM0  (t  - ' c )j 
4 *|  r - _sn<t> | 


(3) 


If  we  now  assume  that  the  observation  point  is  in  the  Fraunhofer  zone  of  the  array 
of  sources,  we  can  approximate  (3)  as 


u(r,  t 


exp 


k (■-»)]  3 

4»r  *-• 


ane 


Ik  • r <t) 
-o  — n 


n=-N 


1 


(4) 


where  r is  the  distance  from  some  reference  point  O (chosen  roughly  at  the  cen- 
o 

ter  of  the  non-vibrating  array)  to  the  observation  point,  as  shown  in  Figure  1.  ^ 
Also _rn  is  the  distance  from  0 to  the  nth  point  source  and  Jjo  = «0/c  k where  k 
is  a unit  vector  in  the  direction  of  the  observer.  We  shall  next  decompose^  as 


6 


• -4 


-n 

Figure  1.  Three  Dimensional  Array  of  Vibrating 
Sources 


r = R + 6 (t) 
— n — n — n 


(5) 


where  R is  the  location  of  the  n**1  point  source  in  the  absence  of  vibration  and 
6n<t)  is  the  perturbation  in  its  position  caused  by  vibrations.  We  can  then  write 


u(r,t)  « Re 


I anexP  tlV  t5n+in(t)1} 
n=-Nj 


(6) 


i* 


I 


We  can  also  obtain  results  for  a continuous  source  distribution.  For  example,  if 
the  unperturbed  antenna  is  planar  we  get 
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JJ  dx  dy  a(x,  y)  exp  {ikQ  • [R(x,  y)  + Six,  y,  t)]  } . 


We  will  now  study  (6)  and  (7)  for  some  special  cases. 


3.  SMALL  HARMONIC  VIBRATIONS 

If  we  assume  that  (<*>0/c)|6n|  « 1 we  can  expand  exp  (ikQ  • 6^)  in  a Taylor 
series.  Also  if  in<t)  is  assumed  to  be  a harmonic  (periodic)  vibration  we  may 
write 


in(t)  = X in | exP  U*  «!*} 


where  w.  is  the  lowest  vibration  harmonic,  6 . = 6*.  because  6 is  real  and 

A n,  -i  -ni  -n 

in0  - 0 because  6n  is  assumed  to  be  the  displacement  about  the  average  position 
Hn-  W we  exrvnd  exp  (ikQ  • 6n)  in  a Taylor  series,  retaining  only  the  first  two 
terms,  substitute  (8)  and  then  take  the  real  part  of  the  result  we  obtain  from  (6) 


<£,  0 = X Unl  c°*  h («-£)  + *S0  • «n  + *n  ] 


2 oo 


• X X |aJ(!V  lin|l)3in  (“o  + i"l)t  + 0nf]  • 


n=-Nj  I =-«o 


where  v(r,  t)  = 4*r  u(r,  t),  <(,  is  the  phase  of  a , 3 , = 0 + k • R + 
o ~ n n n*  n — o — n 

®n t ~ w0r0/c  and  is  phase  of  6^. 
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We  next  assume  that  the  radiation  is  measured  by  the  receiver  system  shown 

in  Figure  2.  For  mathematical  simplicity  we  shall  assume  that  the  filter  passes 

all  frequencies  from  <*>o  - to  “o  + and  completely  rejects  all  others. 

In  this  case  the  output  of  the  filter  is  given  by  (9)  except  with  the  summation  on  f 

running  from  -Mj  to  M2  instead  of  from  -°oto  °o.  Also  if  the  integration  time  T is 

such  that  T » 1 /«■>,,  and  T » l/<o  the  output  of  the  integrator  (ignoring  constants 
'1  ' o 

of  proportionality)  when  the  receiver  is  at  (©,<(>)  is 
T N2  N2 


p <M>  f v2  dt  =}  I I laJ 

o n=-N1  s = -N1 


a cos 
1 s 1 


N2  n2  m2 


(10) 


[*„-*s+V  Rn-k0'5s]  +i  III 

n=-N1  s = -N^  f=-lV 

|anl  I a3 1 (-0  lijMV  I— si  I ) 003  ^ni  'V  • 


The  result  in  (10)  is,  of  course,  the  time-averaged  power  pattern.  This  quantity 
is  more  conveniently  written  as 


P (0,<M  = y 


N2 

I 

n=-N, 


ane 


ik  • R 
— o — n 


M, 


-I  f 
2 ^ 
I =-M 


Nn 


‘ ik  • R 

I a (k  • 6 ) e 

n -o  -ni 

n=-N, 


(11) 
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Figure  2.  Assumed  Receiver  System  (Located  at  0,  <t>) 


9 


The  first  term  in  (11)  is  the  time-average  power  Pq(0,  <t>)  which  would  be  received 
in  the  absence  of  vibration.  Therefore,  the  second  term  in  (11)  represents  the 
correction  to  the  received  power  due  to  the  vibration  of  the  point  sources. 

It  is  interesting  to  examine  (11)  in  some  special  limiting  cases.  Let  us  first 

assume  that  all  the  sources  vibrate  with  the  same  amplitude  and  phase,  so  that 

U I = U.l  and  3>  , = $..  In  this  case  (11)  becomes 
- ni ' '-1 1 ni  l 

p = PD(e,  </>)Ii  + g(e,  <^)1  , (12 


where 


Po<6U>  -i 


N, 


2 

z 

n=-N, 


ane 


ik  • R 
-o  — n 


(13) 


m 2 

g(e,  </>)  = 1 + Z 1^0  ' -ni  I2  • (14 

1 =-Mj 

The  result  in  (12)  means  that  the  locations  of  the  nulls  in  Pq(0,  <t>)  are  not  affected 
by  a uniform  vibration  in  which  all  of  the  sources  vibrate  in  phase. 

A more  interesting  result  occurs  when  the  point  sources  do  not  vibrate  in 
phase,  but  rather  the  phase  of  the  vibration  of  each  point  source  is  random  and 
pncorrelated.  If  the  /th  harmonic  phase,  4>ni,  of  the  vibration  of  the  n source 
#is  independent  of  the  phase,  3>s/,  of  the  sth  source  we  have  that  the  compound 
probability  density  is 


P2(*ni’*si>  = pl(*ni>  p2(*s P 


(15) 


If  we  further  assume  that  the  phases  are  uniformly  distributed,  so  that 
Pji*)  = l/2»  for  0 < $ < 2»,  we  obtain  from  (11)  for  the  ensemble  averaged 
received  power 


<P(0.*>>  = P0<9.*> 


m2 

z 

/ =-Mj 


n = -N 


1 


(16) 
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Finally,  if  in  addition  to  the  random  phase,  the  direction*  of  the  perturbation 
is  also  random  we  must  average  (16)  over  all  directions.  Upon  using  the  result 


2 ti  ir 


j"  d<t>'  j sin  e'  de'  (k  • 6)2  - -y  k^ 


0 0 


where  O’  is  the  angle  between  _k  and  6^  we  find 

k2  M2  N2 

<«».*>>  - Po(0.4)  + -#-  I I laJ2l{ 

i=-Mj  n=-N1 


Xfl2 


From  (18)  it  is  clear  that  if  the  peak  of  the  radiation  pattern  is  assumed  to  be 
at  6 - <t>  ^ 0,  the  average  sidelobe  level  produced  by  the  antenna  vibrations  is 


Z l»„l2  I 


k n = -N! 
o 1 


I 


Therefore,  when  the  phases  are  random  and  uncorrelated,  the  vibrations 
produce  an  average  sidelobe  level  which  can  fill  in  the  nulls  in  the  radiation 
pattern  Pq(0,  <t>),  and  also  alter  the  relative  sidelobe  level  for  the  case  of  a very 
low  sidelobe  antenna.  As  is  clear  from  (19),  the  amount  of  null  filling  will  depend 
on  the  aperture  taper  (that  is,  the  distribution  of  an)  and  the  amplitude  in  wave- 
lengths of  the  vibrations. 


We  assume  that  all  the  sources  vibrate  along  the  same  axis,  but  that  the  align- 
ment of  that  axis  is  completely  random  with  respect  to  the  vector  k . 
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We  can  also  express  (19)  in  terms  of  the  time  averaged  vibration  amplitude, 
provided  the  receiver  bandwidth  is  large*  enough  to  pass  all  of  the  vibration 
harmonics.  Upon  recalling  that 


m / *:<««■  z 

0 I =-« 


where  ^n(t)  is  the  vibration  amplitude  of  the  n radiator,  we  may  rewrite  (19)  as 


i ‘Ski2 


k n=-N, 
o 1 


L ' an 
n = -N, 


As  an  example  of  the  application  of  (21)  suppose  we  have  a linear  array  with 
uniform  amplitude  weighting,  |an|  = 1,  and  the  vibration  is  monochromatic,  so 


6n(t)  = (2)  6RMS  COS  "l  ‘ 


OT  __1  (ko6RMS) 

' 3 Nj  + N2  + 1 


where  + N2  + 1 is  the  total  number  of  elements  in  the  array.  For  example,  a 
21  element  array  with  a 1/20  wavelength  rms  deflection  due  to  vibrations  will  have 


That  is,  provided 


I2,  I I*. 
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a -28  dB  average  sidelobe  level  due  to  the  vibrations.  For  a 21  element  symmetric 

2 ♦ 

(N.  = N0)  array  with  taper  cos  (nir/20)  the  average  sidelobe  level  is  -26  dB. 

1 * tjj 

The  results  we  have  obtained  here  assume  that  the  vibration  of  the  n radia- 
tor is  uncorrelated  with  that  of  the  sth  radiator.  The  analysis  for  the  case  when 
they  are  correlated  is  given  in  Appendix  A. 

4.  LARGE  HARMONIC  VIBRATIONS 

The  results  in  Section  3 are  valid  provided  the  deflections  are  very  much 
smaller  than  a wavelength.  When  this  condition  is  not  satisfied  it  is  not  permis- 
sible to  expand  exp  (ikQ  • 6n)  in  Eq.  (6)  in  a Taylor  series.  In  this  section  we 
will  assume  that  the  magnitude  of  the  deflections  are  arbitrary,  but  we  will  require 
that  the  vibration  be  monochromatic,  so  that 


j>n(t)  = An sin  (tOjt  + 4>n)  , 

where  ^ is  a real  quantity.  In  this  case  (6)  becomes 


(23) 


o 


Nr 


I anexp{iko-  Rn  + ifnsin«o1t  + <tn)} 
n=-N, 


(24) 


where  ?n  = • An>  Upon  using  the  well  known  expansion  of  exp  (ip  sin  <j>)  in 

Bessel  functions,  we  may  rewrite  0 9)  as 


u(r,  t)  = R 


exp 


iu  ft  — — ) * ik  • R ~ im(u,t+$  ) 

^ Y a e ~ ° n T J (€  ) e 1 n 
t-j  n m n 


4^r 


N2 

l 

n=-N 


1 


(25) 


* 2 
This  result  is  obtained  from  (19)  by  using  the  results  in  Jolley  <pp  82-84)  that 


£ cos2  (l&)  : Z cos4  (i§) 

n=-N  n=-N 


3N 
4 * 


2.  Jolley,  L.  (1961)  Summation  of  Series.  Dover,  New  York. 
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If  we  again  assume  that  the  signal  is  passed  through  a filter  which  passes  only 
frequencies  from  to  <*>Q  + ^2^1  ^anc*  rejec^s  all  others),  the  time  aver- 

aged power  pattern  is  given  by 


m2  n2 


P(0,  d>)  = i J u2(r,  t)  dt  =-j-  Y,  Z an  Jm(?n)  eXp  (im<*n  + ‘-o  ' -n) 

tj  m = -M.  n=-N, 


1 “ ”1 


As  was  the  case  before,  if  all  the  sources  vibrate  with  the  same  amplitude 
and  phase  the  radiation  pattern  is  simply  the  product  of  the  nonvibrating  pattern, 
Pq,  and  a correction  pattern.  That  is 


P(8.*>  = I jjL(S)  P(».  <*>)  . 


m - -M , 


where  f = k • A is  a function  of  0 and  <4. 

— — o — 

In  the  other  limit  when  the  phases  $ of  the  vibration  of  each  radiator  are  ran- 

n $ 

dom,  uniformly  distributed  and  statistically  independent,  we  obtain  for  the  ensem- 


ble averaged  power  pattern 


. _ ik  • n 

( P(0,  <t>) ) 4 £ anJ0(?n)e  ° ° 


M2  , N2 


III  • 


m = -M1  n = -Nj 


We  use  the  fact  that  for  n * s 


im(<j  -<6  ) 
<e  > 


2 71 
(A)  f d* n / d*s 


im(d>  -d.  > 
e n s 0 


-Tt  -IT 


In  obtaining  (28),  however,  we  have  not  averaged  over  all  directions  as  we  did  in 
Section  3.  Therefore  (28)  can  be  compared  with  (16)  but  not  (18). 


14 


where  the  prime  denotes  that  the  m = 0 term  is  excluded  from  the  summation. 
When  £n  « 1,  Eq.  (28)  reduces  to  (16),  provided  we  retain  only  the  i = ±1  terms 
in  (16),  so  that  6n  has  only  a single  vibration  frequency. 

If  all  the  radiators  vibrate  in  the  same  direction  and  with  the  same  amplitude, 
so  that  €n  - Cg  = € - ]S0  ' we  cai1  approximate  (28)  by 


m2, 

' N2 

4 

I 4,«> 

Z Uni2 

m = -Mj 

.n=_Nl 

(29) 


Therefore  in  this  case  the  primary  power  pattern  is  altered  by  the  amplitude 
factor  Jq(S)  ^ 1. 

By  using  (29)  and  (13)  we  then  see  that  the  average  vibration-produced  side- 
lobe  level  relative  to  the  peak  power  (at  the  center  of  the  main  beam)  in  the 
absence  of  vibrations  is 


SL  = 


M, 


Z £«> 

m = -M. 


N„ 


Z U I 

t-.  i n i 

n = -Nt 


Nr 


2) 


Z Unl 

vn=-N, 


(30) 


For  (30)  we  see  that  if  ? a ^ is  of  order  unity  the  average  sidelobe  level 
can  be  quite  large.  For  example  if  f = 1,  and  the  array  taper  is  uniform  so  that 
| an|  =1  for  all  n we  get  (provided  Mj,  M2  » 1) 


SL  ~ 


0.416 

n1  + n2  + ! 


(31) 


Therefore,  for  a 21  element  array  (Nj  + N2  + 1 = 21)  we  find  that  the  vibrations 
produced  an  average  sidelobe  level  of  -17  dB.  Consequently,  any  nulls  in  the 
radiation  pattern  would,  on  the  average,  be  filled  in  to  a depth  of  -17  dB  if  ? = 1. 
For  ? - 3 we  find  the  average  sidelobe  level  have  the  even  larger  value  of  -13.  5 dB. 


5.  DISCUSSION 


Which  formula  should  be  used  in  calculating  the  sidelobe  levels  will  depend 
upon  the  physical  conditions.  If  the  vibration  amplitude  is  small  compared  with 
a wavelength  (kQ{  « 1)  and  we  can  actually  estimate  the  lateral  distance  over 
which  the  vibrations  are  correlated  then  clearly  Eqs.  (A6)  or  (A  12)  in  Appendix  A 
are  appropriate.  However,  if  we  are  not  sure  what  the  correlation  distance  will 
be  then  Eq.  (21)  may  be  most  appropriate  for  estimating  the  vibration  induced 
sidelobe  levels.  If  the  vibration  amplitude  is  of  order  or  exceeds  a wavelength 
then  Eq.  (30)  is  most  useful. 


Appendix  A 


Here  we  shall  remove  the  restriction  that  the  amplitude  and  phase  of  each 

radiator  is  statistically  independent  of  those  of  all  other  radiators.  We  shall 

assume,  however,  that  each  radiator  vibrates  along  the  same  axis  in  space  (but 

that  the  orientation  of  that  axis  is  completely  random),  so  that 

k . 6 . = k 6 , cos  0’,  k • 6 , - k 6 , cos  0',  etc.,  where  O'  is  the  angle 
— o — n(  o n t — o —si  os  / 

between  k and  6 In  this  case  (11)  becomes,  after  using  (17)  to  average  over 
— o — ru 

all  vibration  directions 


2 M2 


N2  N2 


<P<8.*)>  = Po<0.*>  111  anaI<6m6st' 


l j n = -Nj  s=-Nj 


— o" (— n"— s^ 


(Al) 


If  the  variables  6 . and  6 . are  assumed  to  be  random  phasors  with  a uniformly 
n*  s"  ^ ^ 

distributed  phase,  it  can  be  shown  that  {6^  t>a()  - P(  (n,  s)  where  the  correlation 
function  p{  (n,  s)  is  real.  Consequently,  we  can  rewrite  (Al)  as 


1.  Beckmann,  P.  (19 GV)  Probability  in  Communication  Engineering,  Harcourt, 
Brace  and  World,  New  York. 
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<P<0, 


k2 

M2 

+ _o 

> 

+ 6 

i a-M 

-R .) 

1 — » 

£ £ anaIp<(n'S> 


.j  = 


(A  2) 


As  an  example  of  the  application  of  (A2)  consider  a symmetric  <Ni  = N2  = N) 
linear  array  of  radiators  lying  along  the  x-axis  with  each  radiator  separated 
from  the  adjacent  one  by  a distance  d.  If  we  further  assume  that  the  deflections 
are  a stationary  random  process,  so  that  p^(n,  s)  = pj(n  - s)  we  can  show  that 
(A 2)  becomes 

k2  M2  N N 

<P<M»  ^ Po<0.*>  +-jf-  Z Z Z anaspi(n"s) 

n=  -N  s = -N 

(A3) 

• exp  (ikQd(n  - s)  sin  6 cos  <t> } 
we  next  assume  that 

pj(n  - s)  = ( Un|  ! 2)  exp 


(n  - s) 


(A  4) 


that  the  array  weighting  is  uniform,  so  that  an  - 1 for  all  n,  and  that  the  mean 
square  deflections  < | 6 nJ  | 2>  are  the  same  for  all  radiators,  so  that 

<Un/!2>  = <U,I2>- 

Then  upon  using  these  assumptions,  substituting  (A4)  into  (A3)  and  using 
(20)  we  can  show  that  if  N » N the  quantity  ( P<0,  <t> ))  can  be  approximated  as 


*We  assume  that  the  filter  bandwidth  in  Figure  2 is  sufficiently  large  that 


00 


I l«Js=  Z l‘„, 


i2  ■*; 


l = -M 


I =-00 
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k2  (62)v/»NN 

< P(0,  0))  ~ P0(fl.<*)  + — 3 exp 


( k 
J . _2 o 


2 ',2  d2  s in2  0 cos2  0 


I ■ 


(A  5) 


where 


i 

<fi2>  =T  f dt<*2<t>> 


Finally,  if  we  assume  that  the  peak  of  the  unperturbed  radiation  pattern  is  at 
0 = <t>  - 0 we  find  that  the  average  relative  3idelobe  is 


SL  ~ 


/ 2N  \ ( k2  N2  d2  sin2  0 cos2  0 ) 

\-WT~i  ) exP  \ - ~^-2 4 j * 


(A6) 


2 1/2 

As  an  example  of  the  application  of  (A 6)  suppose  (6  ) = (2)  ' Spiyjg  cos  “jt, 
d = X/2,  X = wavelength,  = X/20,  2N  + 1 = 21  and  2Nq  = 5.  Then  (A6)  gives 

for  the  average  sidelobe  level  at  0 = 45°,  0 - 0°  the  result  SL  ~ -55  dB.  However, 
for  0 = 45°  and  0 - 90°  the  sidelobe  level  produced  by  the  vibrations  is  -22  dB. 

We  can  also  specialize  (A2)  to  the  case  of  a continuously  distributed  planar 
antenna.  In  that  case  (A2)  becomes 


.2  M2 


P(0,  <*)  = PQ(0,  0)  + -jr-  Y,  ff  dx  dy  JJ  dx'  dy'  a(x,  y)  a'^x',  y') 


f=-Mj  A 


(A  7) 


pf(x,  y;x',  y')  exp  fikQ  • (r  - r')]  , 


where  £ = (x,  y),  _r'  = (x1,  y ' >,  A is  the  surface  area  of  the  antenna  and  a(r)  is  the 
excitation.  We  can  study  (A 7)  for  the  case  when  the  vibrations  are  a stationary 
random  process,  so  that  pj(£,£')  - p#(r  - r').  In  particular  we  shall  assume  that 


pl 


<6/2) 


exp 


(A8) 


If  a(£)  changes  slowly  in  comparison  with  the  variation  in  the  correlation  function 
2 

Pj(r)  and  A » 13  we  may  approximate  (A 7)  by 
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2 ^2 


(P(9,*)>=P0(^)+T  Z <6*2)  ff  dr-2  'a(^) 


l =-M, 


ff dW  I'?* l-'1} 


(A9) 


where  E = x(x  - x')  + y(y  - y').  The  assumption  that  A » /3  allows  us  to  extend 
the  limits  on  the  5 integration  to  infinity.  The  integral  on  ? can  be  performed  by 
writing  d2?  in  cylindrical  coordinates  as  d2?  = I d?  dii  and  kQ  • = kQ  E sin  0 

cos  (<i  - >/,).  Upon  performing  the  integration  on  d2?  we  then  get 


2 M2 


71  K 

(P (0,<t>))  --  P0(e,<t>)  +-g-2-  Z < 6t 2) 


I = -M , 


exp 


^ k2  g2  sin2  e ) 


// 


a(r) ! “ d2r 


(A  10) 


If  we  assume  that  the  filter  passes  nearly  all  the  harmonics  so  that 


Mr 


Z <*/>*  Z <«|2)  =4  f <*2<t»  dt  = <62> 

n r t _ _ * 


I = -M 


I - -no 


(All) 


we  can  use  (A  10)  to  give  for  the  average  sidelobe  level 


SL  ~ ~ k2  (a2)  32 


iJ*  I a(£’ 


a(r)l  2 d2r 


// 


a(r)  | 


exp 


(-5 


k2  fi2  sin2  0 j (A  12) 


For  a uniform  excitation,  a(r)  = 1,  and  0 = 3X,  A * 100*,  0 = 10  and  a single 
harmonic  vibration,  so  that  <a2)  = 26p^lgcos2  <*>jt,  (62)  = ^RjviS’  we  ^'nd  an 
average  sidelobe  level  of  -29.  5 dB. 
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